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Abstract 

Little is understood about how environmental heterogeneity influences the spatial 
dynamics of sexual selection. Within human-dominated systems, habitat modifi- 
cation creates environmental heterogeneity that could influence the adaptive value 
of individual phenotypes. Here, we used the gray catbird to examine if the ecologi- 
cal conditions experienced in the suburban matrix (SM) and embedded suburban 
parks (SP) influence reproductive strategies and the strength of sexual selection. Our 
results show that these habitats varied in a key ecological factor, breeding density. 
Moreover, this ecological factor was closely tied to reproductive strategies such that 
local breeding density predicted the probability that a nest would contain extra-pair 
offspring. Partitioning reproductive variance showed that while within-pair success 
was more important in both habitats, extra-pair success increased the opportunity 
for sexual selection by 39% at higher breeding densities. Body size was a strong 
predictor of relative reproductive success and was under directional selection in 
both habitats. Importantly, our results show that the strength of sexual selection 
did not differ among habitats at the landscape scale but rather that fine-scale vari- 
ation in an ecological factor, breeding density, influenced sexual selection on male 
phenotypes. Here, we document density-dependent sexual selection in a migratory 
bird and hypothesize that coarse-scale environmental heterogeneity, in this case 
generated by anthropogenic habitat modification, changed the fine-scale ecological 
conditions that drove the spatial dynamics of sexual selection. 



Introduction 

The effect of ecological change on evolutionary processes and 
their dynamic interplay remains a fundamental conceptual 
focus in science (Schoener 2011). Moreover, while the relative 
role of ecological factors on evolutionary processes is fairly 
well established for some taxa (e.g., plants), our understand- 
ing of how environmental heterogeneity impacts phenotypic 
selection in most animal populations is limited (Garant et al. 
2007; Kasumovic et al. 2008; Parker et al. 201 1). Quantifying 
these relationships is particularly pertinent given recent evi- 
dence that human activities can alter evolutionary processes 
(Shochat et al. 2006; Smith and Bernatchez 2008). Of particu- 
lar utility will be environment-specific estimates of selection 



that can quantitatively assess the spatial scale of adaptation 
and the mechanistic role of ecological variation on evolu- 
tionary processes (Svensson and Sinervo 2000, 2004). Despite 
the fundamental importance of an ecological-evolutionary 
interplay and evidence that the environmental heterogene- 
ity associated with anthropogenic habitat modification will 
affect the direction and strength of selection (Smith et al. 
2008), few studies have documented how selection varies for 
species that persist in human-dominated landscapes (Smith 
and Bernatchez 2008). 

To date, urban ecology research has primarily focused on 
changes in community structure and the impacts of an- 
thropogenic activities on ecological processes (McDonnell 
and Pickett 1990; Grimm et al. 2008). Moreover, while our 
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understanding of the mechanisms that decouple fundamen- 
tal patterns and evolutionary processes in human-dominated 
landscapes are in their infancy, they have begun receiving 
substantial attention (Shochat et al. 2006). For example, re- 
cent empirical work has documented that the novel selection 
pressures associated with urban environments can result in 
behavioral modifications (Slabbekoorn and Peet 2003), mor- 
phological divergence (Yeh 2004; Smith et al. 2008), diver- 
gent mate choice and reproductive tactics (Candolin et al. 
2007; Perlut et al. 2008), and changes in population genetic 
structure (Bjorklund et al. 2010). Despite these advances no 
study, to our knowledge, has quantified how human habitat 
modifications and the ecological factors associated with these 
changes affect the dynamics of selection on phenotypic traits. 

Here, we focus on sexual selection, which since its incep- 
tion has largely been viewed as a subset of natural selection 
driven by competition over mates and subsequent variance in 
reproductive success (Andersson 1994). Sexual selection can 
operate via female choice (intersexual) and/or male-male 
competition (intrasexual) and is well established as the pri- 
mary mechanism through which conspicuous dimorphism 
in secondary sexual traits arise (Andersson 1994). The op- 
portunity for sexual selection to act is proportional to the 
variance in mating success (Wade and Arnold 1980; Arnold 
and Wade 1984a), and recent evidence suggests that both the 
strength of selection and the expression of secondary sex- 
ual traits can vary both spatially and temporally (Gosden 
and Svensson 2008; Cornwallis and Uller 2010). Tradition- 
ally, male reproductive success was measured as the number 
of social mates a male acquires during the breeding season 
(Webster et al. 1995; Whittingham and Dunn 2005), yet addi- 
tional sources of variation can arise from biased operational 
sex ratios in which not all males obtain partners (Price 1984; 
Dearborn et al. 2001) and from differences in female quality 
such that some females produce more offspring than others 
(Kirkpatrick et al. 1990). A third source of variation stems 
from evidence that social and genetic mating systems are 
often incongruent, especially in birds, with some males sir- 
ing offspring in broods raised by other males (Hasselquist 
and Sherman 2001; Griffith et al. 2002). The offspring re- 
sulting from these extra-pair fertilizations could potentially 
have considerable effect upon variance in mating success and 
subsequent sexual selection, since they increase the success 
of one male while simultaneously decreasing the success of 
another (Westneat et al. 1990). As such, the actual reproduc- 
tive success of males in socially monogamous species is the 
combined total offspring sired from both within-pair and 
extra-pair fertilizations (Webster et al. 1995). 

Given that ecological constraints (i.e., spatial and temporal 
distribution of mates) have long been thought to impose lim- 
its on the degree of sexual selection (Emlen and Oring 1977), 
here we ask how that relationship is modulated by anthro- 
pogenic habitat change. Specifically, we use a two-year study 
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of ecology and paternity in the gray catbird (Dumatella caro- 
linensis) to quantify reproductive strategies and the strength 
of sexual selection on phenotypes important in male-male 
competition and female mate choice in a human-dominated 
landscape. Gray catbirds were ideal for this study because 
they are a socially monogamous species that persists along 
the urbanization gradient yet show substantial variation in 
abundance among habitat types. Moreover, catbirds breed 
in a variety of anthropogenically modified habitats (i.e., the 
suburban matrix [SM] and suburban parks [SP]) that vary 
in habitat quality and breeding density. In particular, density 
has been shown to explain intraspecific variation in extra- 
pair paternity (Gibbs et al. 1990; Gowaty and Bridges 1991; 
Lifjeld et al. 1991; Yezerinac et al. 1999; Stewart et al. 2010) 
and the strength of selection (Bonenfant et al. 2003; Kokko 
and Rankin 2006; Gosden and Svensson 2008). Overall this 
study aims to understand how anthropogenic habitat mod- 
ifications create a cascade of effects on ecological factors, 
reproductive strategies, and the adaptive value of traits under 
selection. 

Methods 

Study species and site 

Our research focused on the gray catbird, a medium-sized 
migratory passerine (ca. 35 g) bird in the family Mimidae 
that breed from Canada to the southeastern United States 
and winters along the Gulf coast to the Caribbean and south 
to Panama (Cimprich and Moore 1995). Throughout their 
range catbirds are found breeding in dense shrubs or vine 
tangles and are often associated with early serai stage suc- 
cessional habitats (Zimmerman 1963; Peck and James 1987). 
These habitat preferences make catbirds common residents 
in both the matrix and suburban parks, yet breeding density 
varies markedly among our study sites with the matrix hav- 
ing lower density than parks (see Results). The catbird mating 
system is traditionally described as social monogamy (Cim- 
prich and Moore 1995) although polygyny has been reported 
in two separate instances (Johnson and Best 1982). Like many 
other monogamous species, catbirds lack distinct plumage 
dimorphism although males tend to average more reddish 
edging on the undertail coverts (Suthers and Suthers 1990). 
Although catbirds lack apparent sexual plumage signals they 
are known to exhibit sexual size dimorphism (Raynor 1979; 
Lent 1990), with males being approximately 3% larger than 
females in several morphological measures. 

We conducted this research from May to August in both 
2008 and 2009 at five study sites (two SM sites and three 
SP sites embedded within the matrix) located in the greater 
metropolitan area of Washington, DC and Newark, Delware. 
Our SM sites were located in Takoma (35.3 ha) and Sil- 
ver Spring, Maryland (18.7 ha), which had human popula- 
tion densities of 1308.9 and 1249.2/km 2 , respectively (United 
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States Census Bureau 2000). Our SP sites were located at 
Sligo Creek Stream Park (9.0 ha), Wheaton Regional Park 
(19.9 ha), and the University of Delaware Ecology Woods 
(16.2 ha) located in Takoma Park and Wheaton, Maryland 
and Newark, Delaware, respectively. Matrix sites were com- 
prised of individual properties within suburban neighbor- 
hoods that varied in size and landscaping but were broadly 
characterized by a broken canopy of large hardwood trees 
(oaks, Quercus; maples, Acer; and sycamores, Platanus) and 
yards that contained a mix of native and nonnative shrubs 
suitable for nesting. In contrast, SPs represented larger con- 
tiguous patches of habitat embedded within the SM charac- 
terized by a closed canopy of hardwoods and a continuous 
understory of dense and largely nonnative shrubs (e.g., mul- 
tiflora rose and Japanese honeysuckle). 

Field methods 

During both years of study, we captured adult catbirds using 
both systematic netting and target netting. Each captured in- 
dividual was banded with a unique combination of colored 
leg bands and an aluminum band from US Fish and Wildlife 
Service (USFWS). At the time of capture a small blood sample 
(20-50 /xL) was taken from the brachial vein and preserved 
in lysis buffer (Longmire et al. 1988) for paternity analy- 
ses. Detailed morphological measurements including wing, 
tail, tarsus (mm), and mass (g) were collected at four of the 
five sites (Delaware excluded) by a single measurer to min- 
imize variation. Nests were located using parental behavior 
and systematic searches of the habitat. We visited each terri- 
tory until we located an active nest and rotated among plots 
to ensure equal sampling effort. Nest contents were mon- 
itored every — two to four days and nestlings were banded 
with USFWS aluminum bands and bled (12-15 /xL) two to 
three days prior to fledging. Despite differences in habitat 
types, nest survival (DSR, daily survival rate) did not differ 
significantly between study plots (/ 2 = 0.52, df = 3, P = 
0.914) or among habitat types (x 2 = 0.05, df = 1, P = 0.814; 
T. B. Ryder, unpubl. data). Social parents of each nest were 
primarily determined using target netting at each nest; how- 
ever, identity of adults was confirmed using observations of 
parental feeding. 

Genetic analyses 

DNA was extracted and purified using a Qiagen BioSprint 
96 robotic system and DNeasy blood and tissue kit (Qia- 
gen Inc., Valencia, CA, USA). We screened 16 microsatellite 
primers developed from Northern mockingbird, Mimuspoly- 
glottos (Hughes and Deloach 1997) and gray catbird (Cabe 
et al, unpubl. data) and optimized the six most polymor- 
phic of these loci for genotyping. Polymerase chain reactions 
(PCRs) were run in 10 /xL volumes consisting of 30 ng of 
genomic DNA, 1 /xL of dNTPs (ImM deoxynucleotide tripo- 
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hosphates), 1 /xL of 5x Taq buffer (2.0 fiM of MgCl 2 , 10 mM 
of Tris-HCl, 50 mM of KC1), 1 fiL of 25 mM MgCl 2 , 0.5 /xL 
of forward and reverse primers, 1 /xL of 2.5 x bovine serum 
albumin (BSA),and0.1 /xL of FlexiTaq polymerase (Promega 
Corp., Madison, WI, USA). PCR products were tagged using 
fluorescently labeled forward primers (Applied Biosystems, 
Inc., Carlsbad, CA, USA). PCR conditions consisted of initial 
denaturing at 94°C for 2 min, followed by 30-40 cycles of 
denaturation at 94°C for 30 sec, X°C annealing for 15-30 sec, 
extension at 72°C for 30 sec, followed by a final extension at 
72°C for 10 min. The annealing temperature (X) was 50°C for 
DCA22 and DCD9, 55°C for DCC24, MP24, and MP25, and 
57°C for DCD22. PCR products were pooled and loaded on 
a ABI 3 1 00x1 automated capillary sequencer together with an 
internal ladder (ROX 500 bp, Applied Biosystems, Inc.). Frag- 
ment sizes and genotypes were assigned using Genemapper 
4.1 (Applied Biosystems, Inc.). 

Each individual was typed at six microsatellite markers that 
varied in the number of independently assorting alleles and 
polymorphic information content (Table 1). We determined 
genotypes for 98% of the total possible cells, then calcu- 
lated allelic frequencies and tested for linkage and departures 
from Hardy- Weinberg equilibrium using Genepop (Rousset 
2008) . None of the loci were linked or out of Hardy- Weinberg 
equilibrium and all loci showed low frequencies of both null 
alleles and typing error (Table 1). The six polymorphic loci 
yielded a high combined probability of exclusion for the pa- 
ternity analyses (P et = 0.99). 

To minimize assigning offspring that match males by 
chance, we used the maximum-likelihood approach avail- 
able in the program CERVUS v. 3.03 (Marshall et al. 1998; 
Kalinowski et al. 2007). Our preliminary simulations used 
10,000 cycles and the true typing error (0.015), measured 
via known mother-offspring comparisons. We used females 
observed at the nest as putative mothers unless they had 
greater than one mismatch in which case the mother was as- 
sumed unknown. Offspring were assigned using both strict 
95% and relaxed 80% confidence as well as a total evidence 
approach (Prodohl et al. 1998; Webster et al. 2004). Un- 
der all three scenarios, paternity was only assigned if the 
presumptive father matched the offspring at five or more 
loci. Using the total evidence approach, we assigned offspring 
only if other nestlings in the same brood were assigned with 
confidence to a sire known to be the attending male at the 
nest. 

Statistical analyses 

We assessed local breeding density (hereafter just breeding 
density) as a potential ecological factor that could influence 
rates of extra-pair paternity and sexual selection. Specifi- 
cally, we used nest locations as a proxy for breeding density 
(pairs/ha). We recorded nest locations (removing renests to 
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Table 1. Details of the six polymorphic microsatellite markers used for gray catbird paternity analyses. 



Locus 


K 1 


N 2 


Hobs 


Hexp 


PIC 3 


Exc 4 


H-W 5 


FNull 


Error 6 


DCA22 


41 


956 


0.955 


0.951 


0.948 


0.819 


ns 


-0.002 


0.012 


DCC24 


30 


953 


0.915 


0.914 


0.907 


0.706 


ns 


-0.001 


0.011 


DCD9 


37 


946 


0.930 


0.946 


0.943 


0.802 


ns 


+0.008 


0.003 


DCD22 


13 


956 


0.871 


0.876 


0.863 


0.599 


ns 


+0.002 


0.033 


MP24 


20 


955 


0.865 


0.885 


0.874 


0.623 


ns 


+0.012 


0.026 


MP25 


28 


890 


0.939 


0.942 


0.938 


0.787 


ns 


+0.001 


0.007 



1 Number of independently assorting alleles. 
2 Number of individuals typed. 
'Polymorphic information content. 
"Exclusion probability of the first parent. 
5 Hardy-Weinberg equilibrium tests. 

6 Error rates calculated from known mother-offspring comparisons. 



avoid overestimating density) for all breeding pairs within 
the study sites using a Garmin etrex Vista HCx and mapped 
those locations in ArcMap V10 (ESRI, Redlands, California). 
To account for fine scale spatial differences in the density of 
breeding pairs, we generated a kernel density surface using 
the spatial analyst extension in ArcMap and the point layer of 
nests as our input. We then generated a continuous measure 
of density to be used in subsequent analyses by extracting 
an individual density estimate for each male from the kernel 
surface. We compared breeding density across sites using an 
analysis of variance (ANOVA) and used binary logistic re- 
gression to determine if density influenced rates of extra-pair 
paternity. 

We calculated the opportunity for sexual selection in the 
two habitat types (SM and SP) as the standardized variance in 
male fertilization success following Webster et al. (1995). The 
standardized variance divided by the square of mean success 
is a measure of the maximum possible strength of selection 
(Arnold and Wade 1984a). We first calculated variance in 
apparent male success (var T A , the number of social young 
produced) and variance in realized reproductive success (var 
T, the number of offspring sired based on molecular pater- 
nity). We second partitioned total variance in realized male 
success into its component parts of within-pair (var(W)) 
and extra-pair fertilizations (var(E)) as well as the covari- 
ance (cov ( W, E)) between the two to determine the effect of 
extra-pair young on the opportunity for sexual selection: var 
( D = var( W) + var (E) + 2 cov( W, E) (Webster et al. 1995). 
Finally, given that comparisons of variance (apparent vs. rela- 
tive) can be misleading about the effect of extra-pair paternity 
on the opportunity for sexual selection (Webster et al. 2007), 
we calculated a Batemans gradient (sexual selection gradient; 
Arnold 1 994) by regression the number of mates against male 
reproductive success. The slope of this relationship provides 
a measure of the intensity of sexual selection arising from 
mating with additional females (extra-pair mates) outside 



the pair bond socially monogamous species (Arnold 1994; 
Arnold and Duval 1994). 

We assessed sexual selection by regressing the number of 
offspring sired against male morphological traits (wing, tail, 
tarsus, mass, and body size) using stepwise least squares re- 
gression. We first used principle component analysis (PCA) to 
collapse the four morphological measurements into a single 
measure of overall body size. The largest portion of variance 
in male body size (45.2%, eigenvalue = 1.81) was accounted 
for by wing and tail, both of which loaded positively on the 
first principle component (PCI). Tarsus and weight loaded 
positively on PC2, which explained an additional 27.3% of 
the variance (eigenvalue = 1.09). Both axes were used as 
predictors in analyses. Second, we calculated linear selection 
differentials (fi) for each habitat using regression coefficients 
(Lande and Arnold 1983). To ensure our selection differen- 
tials were comparable to previous work, we used standardized 
trait values (mean = 0, unit variance) and relative reproduc- 
tive success (number of offspring sired/population mean). 
If individuals were present in more than one year we used 
mean breeding density and cumulative relative fitness over 
both years of study to avoid pseudoreplication in the selection 
analyses. Given the continuous nature of our fitness metric, 
significance of the selection coefficients were assessed using 
the least-square regression analyses. Fitness differentials were 
visualized by fitting relative fitness against traits of interest 
using a cubic spline model and smoothing parameter selected 
via generalized cross-validation (GCV) (Schluter 1988). Cu- 
bic spline models that minimized GCV were fit using a gener- 
alized additive model in program R using the mgcv package 
(R Team 2008; Wickham 2009). The contour surface was es- 
timated using a Loess-smoothing algorithm with polynomial 
regression in Sigmaplot v. 10.0. 

To examine how selection varied by habitat type and an 
ecological covariate, we conducted two analyses using a gen- 
eralized linear model framework. Generalized linear models 
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Table 2. Population-level patterns of extra-pair paternity estimated for gray catbirds nesting in the suburban matrix (SM) and suburban park (SP) sites 
in greater Washington, DC. 





Habitat 


No. of 


No. of 


Percentage broods 


Percentage EPY 


Site 


type 


broods 


offspring 


EPY (W) 


(W) 


Silver Spring 


SM 


22 


62 


22.7 (5) 


9.7 (6) 


Takoma 


SM 


36 


103 


22.2 (8) 


12.6 (13) 


Sligo Creek 


SP 


33 


86 


30.3 (10) 


12.8 (11) 


Wheaton 


SP 


31 


89 


22.6(7) 


13.5 (12) 


UD Ecology 


SP 


43 


115 


25.6(11) 


14.8 (17) 


Total 




165 


455 


24.8 (41) 


13.0 (59) 



used a normal error distribution, and identity link function 
and the maximum likelihood method for estimating model 
parameters. First, we added habitat as a categorical variable 
into the model and examined the interaction between signif- 
icant morphological predictors of male reproductive fitness 
from previous analyses to determine the importance of an- 
thropogenic habitat modification on selection dynamics. We 
specifically asked if the selection differentials observed in 
parks were different from those observed in the matrix by 
comparing the slopes of the two lines. Second, we examined 
both coarse-scale habitat and the fine-scale ecological covari- 
ate by using both habitat context (categorical) and breeding 
density (continuous; pairs/ha) as predictors in the model. 
Again, we were interested in the interaction between density 
and traits under selection. Finally, we compared the mor- 
phology of successful extra-pair sires against the social males 
of those nests in which the extra-pair young were present 
using paired t-tests. All statistics were done using JMP v.8.0.1 
and means and standard errors are reported. 

Results 

Patterns of extra-pair paternity and 
breeding density 

Over two breeding seasons, we sampled a total of 455 off- 
spring from 165 broods at our five sites (Table 2). Of the total 
offspring we assigned paternity to 425 (93.4%), of which 
401 (88%) were assigned at 95% strict confidence, 15 (3.3%) 
were assigned at 80% relaxed confidence, nine (2.0%) were 
assigned using total evidence, and the remaining 30 (6.6%) 
went unassigned. Genetic analyses revealed that 13% of cat- 
bird nestlings were the result of extra-pair fertilizations (all 
95% strict confidence assignments) and that 24.8% of broods 
contained extra-pair young (EPY) (Table 2). We found no 
differences in the rates of extra-pair paternity (EPP) among 
years (2008 = 24.2%, 2009 = 26.4%; y} = 0.53, df = 1, P = 
0.46) or study populations (x 2 = 0.81, df = 4, P = 0.94; 
see Table 2). In addition, we found no significant differences 
in EPP rate by habitat type (x 2 = 0.76, df = 1, P = 0.38) 
despite parks having a slightly higher percentage of broods 



with EPY (26.17% vs. 22.45%). We next examined density as 
a potentially important ecological factor influencing rates of 
extra-pair paternity. The study sites showed significant varia- 
tion in breeding density with parks on average having higher 
density than the matrix (F 3M = 3.85, P = 0.01; Fig. 1). 
Breeding density, as measured by the number of pairs/ha, 
had a significant positive effect on the occurrence of extra- 
pair paternity (x 2 = 4.00, df = 1, P = 0.04, odds ratio = 
0.84). 

Opportunity for sexual selection 

Realized reproductive success of male gray catbirds ranged 
from — zero to nine offspring (mean = 2.03 ±0.15) and did 
not significantly differ by habitat context (f = 0.62, df = 213, 
P = 0.53). Likewise, there were no differences between the 
matrix and parks in the number of within-pair (f = 0.97, 
df = 213, P = 0.33) and extra-pair (t = -1.18, df = 213, 
P = 0.24) offspring males sired. Males that sired extra-pair 
offspring, however, had higher reproductive success in both 
habitats (matrix, t = 3.00, df = 42, P = 0.004; parks, t = 5.03, 



4.0 



£ 3.5 
(0 

a 



■ 
□ 
cn 



2.0 



n = 30 



I Parks 
! Matrix 



n = 17 



n = 24 



n = 27 



Sligo Creek Wheaton Silver Spring Takoma 

Study Site 

Figure 1. Breeding density (pairs/ha) by study site and habitat as mea- 
sured using a kernel density surface. Suburban parks (SP) on average 
tended to have more breeding pairs than the suburban matrix (SM). 
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0 



0 1 2 3 4 5 

Total Number of Mates 

Figure 2. The Bateman's gradient for male gray catbirds, Dumatella 
carolinensis, showing the relationship between the number of genetic 
partners each male had and his reproductive success. Circles are sized 
proportional to sample size (small circles n — 1 , largest circle n = 76). 

df = 78, P < 0.0001). The relationship between extra-pair 
fertilizations and male fitness was further evidenced by a sig- 
nificant relationship between the total number of offspring 
and the total number of partners (F 12 i4 = 460.73, P < 0.0001, 
r 2 = 0.68; Fig. 2). Partitioning the total opportunity for selec- 
tion into its component parts revealed that variance in within - 
pair success (var W 86-91%) was the greatest contributor 
to realized male reproductive success in both environments 
(Table 3). Variance in extra-pair success (var E 7-14%) ex- 
plained a smaller portion of standardized variance, although 
the extent varied by habitat ( Table 3 ) . Covariance explained a 
very small portion of the variance in male reproductive suc- 
cess in both habitats but was positive in the matrix and close 
to zero in park environments (Table 3). While we failed to 
find statistical differences in EPP rates among habitats, parks 
had substantially higher total standardized variance result- 
ing in a 39% increase in the overall opportunity for sexual 
selection. 



Selection differentials 

In a step-wise regression, male body size (PCI) was the best 
predictor of relative fitness with larger bodied males siring 
more offspring than smaller individuals (F 1-151 = 23.43, P < 

0. 0001, r 2 = 0.14). We subsequently examined the selec- 
tion differentials for male body size in each habitat. While 
male body size was under strong directional selection in both 
SP and the SM sites (y6 park = 0.42 ± 0.11, P hH7 = 14.89, 
P = 0.0002, r 2 = 0.15; /J matrix = 0.28 ± 0.09, F L60 = 8.73, 
P = 0.004, r 2 = 0.13; Fig 3), we found no evidence that 
the strength of selection differed between habitats (i.e., slope 
of lines) (x 2 = 0.32, df = 1, P = 0.57; Fig. 3). Moreover, 
a second analysis showed that body size (x 2 = 14.4, df = 

1, P = 0.0001), breeding density (x 2 = 92.3, df = 1, P < 
0.0001), and their interaction (x 2 = 7.6, df = 1, P = 0.0058) 
all influence a male's relative fitness while habitat context 
does not (x 2 = 0.18, df = 1, P = 0.67). These results and 
an examination of a fitness surface indicate that large-bodied 
males in the highest breeding density environments had the 
highest relative fitness regardless of habitat context (Fig. 4). 
Finally, sexual selection for larger-bodied males was fur- 
ther evidenced by significant differences in body size (PCI) 
(Table 4) when making pair- wise comparisons between extra- 
pair sires and the within-pair males they cuckolded. 

Discussion 

Sexual selection is an important force in the evolution of 
male sexual signals (Andersson 1994). Moreover, while sex- 
ual selection studies are common in natural environmental 
contexts, little is known about how anthropogenic activities 
affect reproductive strategies and subsequent selection pres- 
sures when species persist in human-dominated landscapes. 
Here, we document how the reproductive strategies of gray 
catbirds varied within two human-dominated habitats. In 
addition, we measured the relative opportunity for sexual 
selection and estimated selection differentials for traits re- 
lated to male fitness. Our results show that rates of extra-pair 



Table 3. The opportunity for sexual selection varies for gray catbirds breeding in two different ecological contexts. Standardized variances in apparent 
(number of social young produced, var (T A )) and realized reproductive success (number of offspring sired based on molecular paternity, var (7~)) were 
calculated and weighted according to Webster et al. (1995). Components of variance were partitioned to determine the contribution of within-pair 
reproductive success (IV), extra-pair reproductive success (£), and the covariance between the two (W, E) across two breeding seasons. 





Suburban matrix (N = 69) 




Suburban park (W = 146) 




Standardized value 1 


Total percentage 


Standardized value 1 


Total percentage 


var (T A ) 


0.872 




1.27 




var (7") 


0.826 


100.0 


1.15 


100.0 


var (W) 


0.750 


90.8 


0.991 


85.7 


var (£) 


0.059 


7.1 


0.163 


14.2 


cov (W, E) 


0.017 


2.1 


0.001 


0.10 



1 Variance/mean 2 . 
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Figure 3. The relative male mating success of gray catbirds in relation 
to body size (PC1) in the suburban matrix (SM) (blue point and line) and 
suburban parks (SP) (red points and line) in greater Washington, DC. 
Selection differentials (fS) show that body size is under strong directional 
selection in both habitats but the slopes of the lines are not significantly 
different. 



paternity did not differ by habitat (matrix vs. park) but, that 
a finer spatial scale ecological variable, local breeding den- 
sity influenced the probability of gaining extra-pair fertil- 
izations and relative fitness. Partitioning reproductive vari- 
ance showed that while within-pair success was the largest 
contributor to a male's reproductive success in both habitat 
contexts, extra-pair fertilizations explained a higher propor- 
tion of variance in the park habitats. Examining which male 
traits were under selection revealed that larger-bodied males 
have higher relative reproductive success but that the strength 
of selection was density dependent. Specifically, our fitness 
surface shows a peak for large-bodied males at the highest 
breeding densities. These results underscore that estimates of 
selection may be biased if research is only conducted in "high- 
quality/high-density" habitats and that studies should ac- 
count for how selection dynamics are influenced by ecological 
factors. 




rpa "v/, aJ 



Figure 4. Relative fitness in relation to body size (PC1) and breeding 
density (pairs/ha) among male gray catbirds breeding in the suburban 
matrix (SM) and suburban parks (SP) in greater Washington, DC. The 
contour surface was estimated using a Loess-smoothing algorithm with 
polynomial regression. 

Patterns of extra-pair paternity and 
breeding density 

Since the discovery that a substantial number of socially 
monogamous species exhibit sexual promiscuity (extra-pair 
paternity), biologists have put forward adaptive (Westneat 
and Sherman 1997; Petrie and Kempenaers 1998; Griffith 
et al. 2002) and nonadaptive hypotheses (Arnqvist and Kirk- 
patrick 2005; Forstmeier et al. 2011) to explain intra- and 
interspecific variation in EPP rates. To date, the majority of 
work at both taxonomic levels has focused on ecological fac- 
tors that can impact the spatial and temporal availability of 
mates yet few studies have examined how these processes 
are modulated by anthropogenic habitat change (i.e., frag- 
mentation and urbanization; but see Kasumovic et al. 2008; 



Table 4. Comparison of male gray catbird morphological traits of extra-pair sires and the males they cuckolded. 







Extra-pair sires 


Cuckolded males 




Trait 


N 


(mean ± SE) 


(mean ± SE) 


P 1 


Wing (mm) 


23 


90.57 ± 0.51 


90.48 ±0.47 


0.90 


Tail (mm) 


22 


95.14 ±0.76 


94.45 ±0.63 


0.49 


Tarsus (mm) 


17 


27.93 ± 0.23 


27.31 ±0.18 


0.04 


Mass (g) 


22 


35.03 ± 0.28 


34.39 ±0.32 


0.14 


Body size (PC1) 2 


17 


0.63 ± 0.21 


-0.04 ±0.1 7 


0.02 



'P-values from paired t-tests. 

2 Wing and tail accounted for the largest portion of variance in male body size. 
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Townsend et al. 2011). Overall, our results show that gray 
catbird populations breeding in the matrix and parks have 
moderate levels of extra-pair paternity (25% of broods, 13% 
of offspring). While we failed to find statistical differences 
in EPP rates among populations or habitat type, breeding 
density did influence the probability of a brood containing 
extra-pair offspring in both habitat types. These results sug- 
gest that fine-scale ecological differences within habitat type 
can affect reproductive opportunities. 

These results are also consistent with comparative evidence 
that density can play an important role in explaining EPP 
variation at the intraspecific level (Westneat and Sherman 
1997; Moller and Ninni 1998; Griffith et al. 2002). Although 
we observed a relationship between density and EPP in our 
catbird populations, the relationship need not be causal if 
other factors that impact EPP rates and breeding density co- 
vary. For example, breeding density may be related to habitat 
structure such that male mate guarding is less efficient in 
the high shrub density environments of park habitats (West- 
neat and Sherman 1997). Higher breeding density may also 
reflect habitat quality and food availability such that males 
breeding in those environments are "emancipated" from a 
certain degree of parental care and can more readily seek 
extra-pair copulations (see Hoi-Leitner et al. 1999). Alterna- 
tively, extra-pair behavior maybe constrained in urban land- 
scapes because of the costs associated with seeking extra-pair 
mates (i.e., greater distances between territories; Norris and 
Stutchbury 2001). Regardless of the exact mechanism, indi- 
viduals breeding in higher density environments experience 
higher mate encounter rates and are more likely to engage in 
extra-pair behaviors. 



Opportunity for sexual selection 

In socially monogamous taxa, differential reproductive vari- 
ance is largely driven by rates of extra-pair paternity, which 
can vary among populations based on ecological or behav- 
ioral factors (Griffith et al. 2002). Given that the strength 
of sexual selection is thought by some to be proportional 
to variance in mating success (Wade and Arnold 1980; 
Arnold and Wade 1984b; but see Klug et al. 2010a), extra- 
pair fertilizations, should theoretically increase the oppor- 
tunity for selection (reviewed in Whittingham and Dunn 
2004). Moreover, while there is little doubt that EPP can 
increase the opportunity for sexual selection, the magni- 
tude of its effect appears to vary significantly among species 
(Freeman-Gallant et al. 2005; Whittingham and Dunn 2005; 
Albrecht et al. 2007; Dolan et al. 2007). In particular, esti- 
mating the effect of extra-pair success on the opportunity 
for sexual selection may be sensitive to population sampling 
(i.e., the proportion of sires sampled; Freeman-Gallant et 
al. 2005). Our cumulative assignment rates for all offspring 
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(93%) and the extra-pair subset (67%) suggest that we ad- 
equately sampled our study populations and that our es- 
timates of the opportunity for sexual selection are largely 
unbiased. 

A second possible challenge associated with understand- 
ing the effects of extra-pair paternity on the opportunity 
for sexual selection has been the reliance on comparisons of 
standardized variances (apparent vs. realized). Recent work 
suggests that variance comparisons can be misleading and 
that extra-pair paternity can contribute to the opportunity 
for selection even in the absence of differences between ap- 
parent and realized success (see Webster et al. 2007; Klug 
et al. 2010a; Krakauer et al. 2011; While et al. 2011). Here, 
our results show a strong positive relationship between the 
number of partners a male has and his reproductive success 
(Bateman's gradient; Fig. 2) despite negligible differences in 
apparent and realized variance. Thus, males increased their 
fitness by engaging in extra-pair behavior and mating with 
females outside the social pair bond. Ultimately, if males 
in a socially monogamous system can gain via extra-pair 
fertilizations and increase their fitness, then the traits as- 
sociated with extra-pair success will be favored by sexual 
selection (Webster et al. 2007). In this way, the additional 
variance associated with extra-pair mating can contribute to 
the strength of sexual selection on specific male phenotypic 
traits. 

To date, most studies of socially monogamous passerines 
have shown that standardized variance in male reproductive 
success is largely explained by within-pair success (Webster 
et al. 1995, 2001; Freeman-Gallant 2005; Whittingham and 
Dunn 2005; but see Dolan et al. 2007). Our results largely 
corroborate this finding in which within-pair mating success 
explains the majority of the standardized variance (86-91%) 
in both habitats. In contrast, extra-pair success explained a 
smaller portion of the variance (7-14%) although this varied 
by habitat type. In the higher breeding density park environ- 
ment, extra-pair success explained considerably more of the 
standardized variance in male reproductive success and was, 
in large part, responsible for the subsequent 39% increase 
in the opportunity for sexual selection. While covariance ex- 
plained a marginal portion of the variance, the positive values 
in the matrix suggest that males were able to gain extra- 
pair fertilization without the trade-off of losing offspring in 
their social nest. Such trade-offs (extra-pair gains = within- 
pair losses) are often used to explain the limited impact of 
extra-pair mating on the opportunity for sexual selection 
(Webster et al. 1995). In contrast to the matrix, covariance 
in the higher density park environment was close to zero. 
These results suggest that males in environments with higher 
male-male competition cannot concurrently maximize both 
within- and extra-pair success because of the potential trade- 
offs between mate guarding and engaging in extra-pair 
behavior. 
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Sexual selection, habitat, and local breeding 
density 

Sexual selection can operate on a phenotypic trait via female 
choice (intersexual) and or male-male competition (intra- 
sexual) and is supported when a trait covaries with repro- 
ductive success (Price 1984; Andersson 1994). While traits 
under sexual selection vary widely among taxa, a number of 
studies have shown body size to be a target of sexual selection 
(Andersson 1994; Gontard-Danek and Moller 1999). Here, 
we report that body size was a predictor of a male's relative 
reproductive success (WPY and EPY) for catbirds breeding 
in both matrix and park sites. While male body size is clearly 
related to reproductive output in catbirds, we are currently 
unable to differentiate the relative roles of intra- versus inter- 
sexual selection on trait elaboration. On one hand, large body 
size could provide advantages during male-male competition 
contests as has been shown for red-winged blackbirds (Searcy 
1979; Eckert and Weatherhead 1987). Alternatively, female 
preference for male size could evolve via direct phenotypic 
benefits in which larger males provide material advantages 
(e.g., territory quality and/or parental care) or via the indi- 
rect genetic coupling of preference and trait (e.g., Fisherian 
and indicator mechanisms; Andersson 1994). Regardless of 
the mechanism, body size is heritable in birds (Garnett 1981; 
Wiggins 1989) such that selection for this male phenotype will 
yield larger offspring that may have a survival advantage in 
both environmental contexts. Ultimately, if large male body 
size affects viability attributes such as survival and reproduc- 
tive success it may be favored by both natural and sexual 
selection in catbird populations. 

While sexual selection may be the most parsimonious ex- 
planation for our findings, an alternative hypothesis is that the 
observed relationship between body size and fitness is indi- 
rectly caused by variation in rearing environment (see Bolund 
et al. 2011). As mentioned earlier, if breeding density varies 
with habitat quality than differences in rearing environments 
could indirectly influence body size even in the absence of 
strong selection. This idea is indeed plausible, yet in addition 
to the strong correlation between relative fitness and body 
size, our results also show that body size influenced extra- 
pair mating opportunities. Specifically, extra-pair males were 
on average larger than the within-pair males they cuckolded. 
This second line of evidence is also consistent with hypothesis 
that sexual selection is acting on male phenotype. Ultimately, 
only methods that can decouple the relative roles of genetic 
and environmental variation on body size (e.g., animalmodel 
approaches) will determine if the relationship between body 
size and fitness is causal. 

Here, we aimed to understand how habitat context and its 
ecological covariates influence the strength of sexual selection 
in gray catbirds. At our study sites, habitat and breeding den- 
sity were largely confounded, with parks having high density 



and the matrix having low density. While both habitat types 
occur within human-dominated landscapes, all of our study 
sites fell well within the range of expected variation in breed- 
ing density (e.g., 0.3 to >10 pairs/ha; Harcus 1973). Our re- 
sults suggest that habitat type does not influence the strength 
of selection, yet the limited overlap in breeding density be- 
tween park and matrix environments may have influenced 
our ability to detect the effect of coarse-scale environmental 
differences. As suggested by Wade and Kalisz (1990) ma- 
nipulations of the ecological agents that drive selection are 
required to fully understand the selection process. As such, to 
fully understand the relative roles of habitat and an ecological 
covariate, such as density, on the strength of sexual selection, 
an experimental manipulation would be required. 

While it has long been established that breeding density 
(i.e., the spatial and temporal distribution of mates) can be 
an important force shaping animal mating systems (Emlen 
and Oring 1977), only recently has work begun addressing 
how density influences the spatial and temporal dynamics 
of selection (Kokko and Rankin 2006; Gosden and Svens- 
son 2008). Thus far, results are mixed ranging from density 
having a strong impact on the strength of sexual selection 
(Conner 1989; Clutton-Brock et al. 1997; Bonenfant et al. 
2003; Kasumovic et al. 2008; Lehtonen and Lindstrom 2008) 
to little or no effect (Head et al. 2008; Klug et al. 2010b). 
Our results show that male body size in gray catbirds was 
under the strongest selection in the highest breeding density 
environments. In contrast, lower density environments likely 
experienced relaxed sexual selection because of limited mate 
choice opportunities and mate acquisition (Candolin et al. 
2007; Candolin and Heuschele 2008). Other recent empirical 
work has hypothesized that localized density- and frequency- 
dependent interactions between the sexes may emerge as a 
primary driver of spatial variation in sexual selection (Gos- 
den and Svensson 2008). Viewed cumulatively our results 
support this hypothesis and suggest higher density environ- 
ments increase the probability of extra-pair fertilizations, the 
opportunity for and the strength of sexual selection. More 
broadly, our results suggest that coarse-scale environmen- 
tal heterogeneity, in this case generated by anthropogenic 
habitat modification, can influence the fine-scale ecological 
conditions that drive reproductive strategies and the spatial 
dynamics of sexual selection. 
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